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Abstract. Dynamics of the Brans-Dicke theory with a scalar field potential function is 

investigated. We show that the system with a barotropic matter content can be reduced 



to an autonomous three dimensional dynamical system. For an arbitrary potential 
function we found the values of the Brans-Dicke parameter for which a global attractor 
2^ ! in the phase space representing deSitter state exists. Using linearised solutions in the 

vicinity of this critical point we show that the evolution of the universe mimics the 
ACDM model. From the recent Planck satellite data, we obtain constraint on the 
variability of the gravitational coupling constant as well as the lower limit of the mass 



5-h of the Brans-Dicke scalar field at the deSitter state. 
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1 Introduction 

Recent astronomical observations [1] indicate that the Standard Cosmological Model 
(ACDM) is still a good effective theory of the physical Universe, although some anoma- 
lies in the power spectrum of CMB are predicted [2]. 

In this theory the crucial role play cosmological parameters which values are 
estimated from the astronomical observations. The nature of some parameters, like 
density parameters for dark sector of the Universe (dark matter and dark energy), is 
still unknown up to today and we are looking for an explanation. We believe that more 
fundamental theory of gravity will reveal the nature of the cosmological parameters. 
In other words we expect that the nature of these parameters is emergent. 

The Brans-Dicke theory of gravity [3] seems to be an interesting way toward 
construction of some more adequate description of evolution of the Universe by which 
one can explain why the universe enters the accelerating phase of expansion just in the 
current epoch. 

In this paper, we show how dynamical evolution of the standard cosmological 
model emerges from the cosmological model in which, instead of General Relativity, 
we use the Brans-Dicke theory of gravity. In principle, we obtain two different evolu- 
tional scenarios in which the ACDM model with some additional corrections is recov- 
ered. Those new additive terms appear in the equation describing the evolution of the 
Universe as well as the parameters describing the dark sector are modified. The gen- 
eralised Friedmann-Robertson- Walker equation can be tested trough the astronomical 
observations and the model itself can be falsified be the data. 

The gravitational sector of the Brans-Dicke theory is modified by the presence 
of some potential function, while the barotropic matter is used for description of the 
matter content of the universe. However, because of the indetermination of the form 
of the potential (for some new results concerning the scalar field potential function for 
the inflaton field see recent Planck results [4]) our strategy is to obtain results without 
the "potential bias" as soon as possible. Therefore in our analysis we concentrate on 
the behaviour of the dynamical system describing the evolution of the universe in the 
vicinity of a special critical point whose position in the phase space does not depend 
on the detailed form of the potential function. 



We demonstrate that there are two different types of the emergence of the ACDM 
model. Different initial conditions and the model parameters give rise to different 
scenarios. In the first scenario the ACDM model is approached monotonically while 
in the second approach we obtain oscillatory approach to the deSitter state. In both 
scenarios, the leading term represents evolution in the ACDM model. 

Only the astronomical data can help in selecting the type of evolution which is 
realised by our universe. Alternatively, deeper knowledge of the initial conditions and 
the model parameters, can distinguish a single type of the behaviour. 

2 The Brans-Dicke cosmology 

The action for the Brans-Dicke theory [3] in the so-called Jordan frame is in the 
following form [5, 6] 

S= /"d 4 x v ^J0 J R-^V a 0V a 0- 2V{4>) \ + lQnS m (2.1) 

where the barotropic matter is described by 

S m = [d 4 x^£ m , (2.2) 



and w BD is a dimensionless parameter of the theory. 

Variation of the total action (2.1) with respect to the metric tensor SS/Sg^ = 
gives the field equations for the theory 



Wbd 



V M V v <j> - l - g^ V a <p V a <p\ - g^V((f>) - (g^Of) - V M V„< 



kT (m) 



(2.3) 



where the energy momentum tensor for the matter content is 



T - , = -^^ ( ^ £m) ' (2 ' 4) 



Taking the trace of (2.3) one obtains 



i?=^V^V Q + 4^ + 3^-^T- (2.5) 



Variation of the action with respect to gives 



□0 = -L V> V Q --$-(R- 2H0)) , (2.6) 

10 ^Wrfi 



and using (2.5) to eliminate the Ricci scalar R, one obtains 

The field equations (2.3) are adopted to the spatially flat Friedmann-Robertson- 
Walker metric 

ds 2 = -d£ 2 + a 2 (t) (dx 2 + dy 2 + dz 2 ) (2.8) 

and the matter content is described by the barotropic equation of state p m = w m p m , 
where p m and p m are the pressure and the energy density of the matter. 
The energy conservation condition is 

oj_bu<P 2 , V{(j>) OTJ <P , 8tt 
2 

and the acceleration equation is 



3H 2 = ^ + ^l-3H^ + ^ Pm (2.9) 



h _ ^ BD 2 1 2V{<f>)-<t>V'{<f>) ^ 2H <f> M p 2 + oo BO (l + w m ) 

2 (f> 2 3 + 2u BD (p <fi <p " 3 + 2u BD 

The dynamical equation (2.7) for the scalar field reduces to 

+ 3#0 = 2 — — + 87rp m — — . (2.11) 

In what follows we introduce the following energy phase space variables 



Then the energy conservation condition (2.9) can be presented as 

I^ + '-T*-*' (2 ' 13 > 

and the acceleration equation (2.10) 

H n __ UJ BD 2 



— =2x- ^x 2 y 2 (2 + A) 

H 2 2 3 + 2u BU y V ; 

„ (. W BD 2 2 \ 2 + W B d(1 + W m ) 

— 3 1 + x — x — y ' 



(2.14) 

-Wm) 

6 'J 3 + 2w BD 

The dynamical system describing the evolution of the Brans-Dicke theory of grav- 
ity with a scalar field potential and the barotropic matter has the following form 

dx o 2 H 6 2 / W BD 2 2^ l-3tUr 

-3x — x — x—— + - — j/ (2 + A) + 3 1 1 + x — x — y 



dr H 2 3 + 2w Rn y 6 * ' / 3 + 2u 

|^-y(i,(i + A) + |], (2,5) 

^ = xA(l-A(r-l)), 



where 4- = -tt— and 
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where ()' = ^. 

From now on we will assume that V = T(A). The critical points of the system 
(2.15) depend on the explicit form of the T(A) function. One can notice that the critical 
points (x* = , y* = ±1 , A* = — 2) do not depend on the assumed T(A). Additionally, 
the acceleration (2.14) calculated at these points vanishes, giving rise to the deSitter 
evolution. 

Our analysis we concentrate on the critical point (x* — 0, y* — 1, A* = — 2) which 
corresponds to the deSitter state, while the critical point with y* = — 1 corresponds to 
the anti-deSitter state. 

The linearisation matrix calculated at this point is 



o 2+2iL) B - D +3Wm —f\ — 



3ui„ 



3+2w BD 3+2w BD 3+2o; B d 

A = I 3 l+2ui BD w m ^ 2+tJBD(l+w m ) 3 I (o 1 7) 

I 2 3+2w B D 3+2o; B d 3+2w B d ' \ ■ ' ) 

|<5 



where we assumed ^~ 

and ^8r 



= A*(l — A*(r(A*) — 1)] = |<5, and 5 is a nonzero constant 
is finite. The eigenvalues are 




l 1 = -3(l + Wm ), l 2<3 = - - 1 ± J BD . (2.18 



For 5 = one of the eigenvalues vanishes and this critical point is degenerated. There- 
fore the analysis based on the linearisation matrix is not enough and the center manifold 
theorem must be used [7-9]. 

The critical point under consideration is stable when w m > — 1 and 

< (2.19) 



3 + 2tu 



for — 1 < ttt < is a stable node, while for ^-n < —1 is a stable focus. This 

3+2o) B D ' 3+2w B d 

indicates that for a given family of potential functions described by the function T(A) 
we can always find ranges of the cu BD parameter for which the critical point under 
considerations is a stable one. Let us consider two simple potential functions. For 
potential of the type V(<f>) = Vo(0 2 — v 2 ) 2 , we have 5 = — y, and for u BD > I we have a 
stable node critical point, while for — | < co BD < | we have a stable focus. It represents 
a saddle type critical point otherwise. For the potential type V((p) = \m 2 (f) 2 + |0 4 , 
we have S = |. For u BD < — ^ this critical point represents a stable node, and for 
— y < u BD < — | we have the critical point of a stable focus type. 
We need to consider both types of behaviour separately. 



In the first case we use the following substitution 



S 4 

-n(n-3) (2.20) 
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where for < n < | we have a stable node critical point. The eigenvalues of the 
linearisation matrix are 

k = -3(l + w m ), l 2 = -n, l 3 = -3 + n, (2.21) 

and the linearised solutions are 

= 4 "("-3)('+^y-^) _ + 

d(n - 3(1 + w m ))(n + 3w m ) 

3o(2n — 3)(n — 3(1 + w m J) \ 

+ (n-3(l + w m ))(35Ax-8(n-3)AA) j exp (-nr) + 

n — 3 / 

+ M (2,^3)(, i + 3u. m ) ( 4 " ( " " 3)(1 - 3t "'" )(AX - 2AS) + 

+ (n + 3w m )(35Ax + 8nAA) J exp ((-3 + n)r) , (2.22a) 

( \ i , ( ~n(n - 3)(1 + w m )(l - 3w m ) IV * , . 

+ eS(2n - 3) (n "- 3(1 + mm) ) ( - 4 "'" - 3 >(' - 3l """ A * - 2Ay » + 



(n — 3(1 + w 



m ))(35Ax - 8(n - 3)AA) J exp ( - nr) + 



n — 3 / 

6o(2n — 3){n + 3u> m ) \ 

+(n + 3u> m )(3<5Ax + 8nAA) j exp ((-3 + n)r) , (2.22b) 

w n n(n — 3)(1 — 3w m ) ,. . . , , N 

a(t) = - 2 - 2 (n -3(i + !L))( n+ L) (Az - 2As) exp ( - 3(1 + Wm)T) + 

4n(n - 3)(1 - 3m m )(Ax - 2 Ay) - 



i(2n - 3)(n - 3(1 + w m )) 



(n - 3(1 + w m ))(35Ax - 8(n - 3)AA) J exp ( - nr) 
' 4n(n - 3)(1 - 3w m )(Ax - 2Ay) + 



i(2n - 3)(n + 3w r , 

+ (n + 3w m )(35Ax + 8nAA) j exp ((-3 + n)r) (2.22c) 



where Ax = x^\ Ay = y^ l > — 1 and AA = A^ + 2 are the initial conditions. 

Using these linearised solution and the conservation condition (2.13) we can find 
linearised solution for the barotropic matter density fraction. Up to linear terms in 
initial conditions we obtain 

Q m « (Ax - 2Ay) exp ( - 3(1 + w m )r) = fi mji exp ( - 3(1 + w m )r). (2.23) 

Now, using the linearised solutions, we are ready to compute the Hubble's func- 
tion. The acceleration equation (2.14) can be cast into the following form 

din if 2 H 2 6 2 / u BD 2 2 \ 2 + cj bd (1 + w. 



2— = 4 x - Ubo x 2 --— y 2 (2+A)-6 (l + x- ^x 2 - y 



dr H 2 au 3 + 2w BD + V 6 "J 3 + 2w BD 

(2.24) 

after integration of this equation, and up to the linear terms in initial conditions we 

obtain 

tt I \ \ 2 / \ — 3(l+w m ) / \ — n / \ — 3+n 



- n Afi + n Mfi — + n ni0 — + n 3nfi — , (2.25) 

H(clq)J \a J \a J \a J 

where 

A 4n(n-3)(l-3^ m .)(4 + 3 Wm ) ^ 

^m,o = 1 - -7777 — -= 77 777- 77" Mm,0 , (2.26a) 

V 3d(n + 3w m )(n - 3(1 + w m )) J 

w + 1 ( ' ln(n - 3)(1 - 3?vj m ) QU A/ao N 

^n,0 = 77777 57 577- v ^m,i ~ 3d Ax + 8fl - 3 AA I — - 

3o(2n — 3) \ n — 3(1 + w m ) J \a^h 

(2.26b) 

n-4 / 4n(n-3)(l-3w m ) n *A\\( a °Y 3+n 

^3n,0 =77777 77 7^ ^m,i ~ 35 Ax - 8nAX ^"7 

3o(2n — 3) \ n + 3w m J Ka^' / 

(2.26c) 



and 

^A,0 — 1 — ^M,0 ~~ ' ^n,0 — ^3n,0 j (2.27) 

and we have used that in linear approximation 

/ a® \ -3(i-H«m) 

^m,0 = ^m,i ( - 77" J • (2.28) 

Now we proceed to investigation of the critical point of the stable focus type. In 
this case we use the substitution 

-(9 + An 2 ) (2.29) 
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and the eigenvalues of the linearisation matrix are 

3 3 

h = -3(1 + w m ) , l 2 = --- ±n, l 3 = -- + ±n. (2.30) 



The linearised solutions of the system are 

,x 4(4n 2 + 9)(l + w m )(l-3m m ) /A , , , 

<r) = 6{4n2 + 9{l + 2Wm)2) (A* - 2 A,) exp ( - 3(1 + w m )r) + 

/ _ 4(4n 2 + 9) (1 + ^) (1 - 3u, w ) _ A _ 3 

V 5(4n 2 + 9(1 + 2w m ) 2 ) v y; / PV 2 y v ; 

1 / 2(4n 2 + 9)(4n 2 - 9(1 + 2w m ))(l - 3w m ) ( _ , 

6n\ 5(4n 2 + 9(1 + 2w m ) 2 ) [ % V '~ 

4 \ 3 

-9Ax - -(4n 2 + 9)AA exp ( - -r) sin (nr) , (2.31a) 

o J 2 

*> - ^ ( 4(4 td^^^r ' - »** - **> - < - 3<i + -^ + 

4 ( - 4(4 t ( ^^^r !) ^ - "»> + **) - < - r> - w + 

1 / 2(4n 2 + 9)(4n 2 -9(l + 2 Wm ))(l-3u> m ) , _ , 
12nV 5(4n 2 + 9(1 + 2u> m ) 2 ) ^ X yj ~ 

-9Ax - -(4n 2 + 9)AA ] exp ( - -r) sin (nr) , (2.31b) 

o ) 2 

^ - - 2 - Srfr^) ^ - 2A *> exp < - 3(1 + •*■") + 

/ (4n 2 + 9)(l-3w m ) , A A , \ , 3 , , , 

+ ( 2(4n 2 + 9(1 + 2 Wm T 2 ) ^ ~ 2A ^ + AA ) 6XP ( - 2 r ) C ° S W + 

3 / _ 2(4n 2 + 9)(l + 2 Wm )(l-3 Wm ) _ 

+ 8nV 4n 2 + 9(l + 2 Wm )2 l Ax z/x y) + 

+5 Ax + 4AA j exp ( r) sin (nr) . (2.31c) 

where Ax = x®, Ay = y® — 1 and AA = A^ + 2 are the initial conditions. 

Using these solutions and the acceleration equation (2.24) one is able to obtain 
the Hubble's function in the vicinity of the critical point of the focus type. Up to linear 
terms in initial conditions we have 



H(a)\ 2 „ „ / a \^(i+ Wm ) 



H(a )J \a 



^A,0 + &M,0 — + 



— ) Qcos,o cos I n In ( — ) ) + Q sin ,o sin ( n In ( — 

a J \ \ \a JJ V V a o 



(2.32) 



where 



^A,0 — 1 — ^A/,0 ~~ &<cos,0 (2.33) 



and 



4(4n 2 + 9)(l-3 Wm )(4 + 3 Wm/ 
Qm '°- X 35(4n 2 + 9(1 + 2^) ^ m '°' (2l34a) 



Q 1 / 4(4n 2 + 9)(l-3u> m )(4 + 3tQ n 

cos '° 35 V 4^ + 9(1 + 2^)2 

-3,A, + 8Aa)(«)^cc.( B 1»(*)). 

1 / 2(4n 2 + 9)(1 - 3^ m )(4n 2 - 15(1 + 2w m )) 
+ Q5n{ 4n* + 9(l + 2w m y m ' J+ 



155 Ax + 4(4n 2 + 15) A\\ f^\ *'* sin (n In f^\ ^ , (2.34b) 



n 1 / 2(4n 2 + 9)(1 - 3w m )(4n 2 - 15(1 + 2u> m )) n 



6Sn 



V 4n 2 + 9(1 + 2w m ) 2 



+ 155Ax + 4(4n 2 + 15)AAJ (-^-) 2 C os(nln 



/ VaW/ """ V'" ~* V a(i) 



1 / 4(4n 2 + 9)(l-3 Wm )(4 + 3tQ n 
35 V 4^ + 9(1 + 2^)2 

-35Ax + 8AA)(-^)- 3/2 S1 n(,ln(^)). (2.34c) 

In this section we obtained two forms of the Hubble's functions for two different 
types of behaviour in the vicinity of the critical point representing the deSitter state. 

2.1 Dust matter 

For dust matter and special initial conditions: 

4 1 

LAX = -ri'rn.i j ^^ ~ T:^^m,,i \Z.OD) 

I 

both forms of the Hubble's functions (2.25) and (2.32) take the same form. Namely, 
in the case of the monotonic approach to the deSitter state (2.25) we have Q n fi = 
and n 3n! o = 0, and in the case of the oscillatory approach (2.32) we have f2 cos0 = 
and fl S in,o = and the resulting form of the Hubble's function is 

|^r V « 1 - ^m,o + ^M,o (-) (2.36) 

H{a Q )J \a J 

where 

^m,o = ( 1 — t^t ] ^m,0 • (2.37) 

If we assume that in the model we include only the baryonic matter f2 m0 = ^bm,o> then 
from the recent astronomical observations of the Planck satellite [1] we have that the 
present total matter density parameter is Qm,o ~ 0.315 and that the baryonic matter 



density parameter is £lbm,o ~ 0.049. This gives us an opportunity to find the value of 
the 5 parameter of the model, which gives us information about the second derivative 
of the potential function at the critical point. After a little algebra one finds 

16 
S = —, r- ps -0.9825 . (2.38) 

3 [!-£*' 



Next, from the definition of the S parameter one can directly calculate the value of the 
second derivative of the scalar field potential at the critical point 



V"{(j>)V{<l>) 



vw 



0.5921, (2.39) 



which indicates that at the critical point the second derivative of the potential function 
is positive V"(0)| > 0. 

Then from the assumed initial conditions (2.35) one gets the present values of the 
phase space variables 

x (a ) == -Jj _ « -0.1995 , A(o ) = -<P-rrnl « -2.0245 . (2.40) 



H<p 



o 



In the Brans-Dicke theory the field <fi is identified as the inverse of the gravitational 

coupling which now gives 

G 

0.1995 (2.41) 



HG 



ii 



which, taking the present age of the universe to = 13.817 x 10 yr, gives 



G 

G 



1.44 xlO" 11 —, (2.42) 

yr 



and is in good agreement with other observational constrains on the variability of the 
gravitational coupling constant [10]. 

One can also calculate the mass of the Brans-Dicke scalar field. In the Jordan 
frame it is given by [11] 



m 



2 



■{WM-V'M), (2.43) 



3 + 2w BD 
which in the previously introduced phase space variables becomes 

6 



m. 



2 



3 + 2u B] 



-H 2 y 2 X(l + XT(X)). (2.44) 



Direct calculation of the asymptotic value of the scalar field mass at the critical point 
gives 

m 2i =J_ul i (2.45) 

'* 4 *3 + 2w Rn v ; 



where H% pa Hq(1 — Qm,o) is the asymptotic value of the Hubble's function and H its 
present value. Inserting the estimated value of the S parameter and Hq fh 1.5 x 10 _33 eV 

one obtains 

10" 33 

m\ sa 1 S4 — e V , (2.46) 

V3 + 2w BD 

which is an asymptotic value of the mass of the Brans-Dicke scalar field at the deSitter 
state. Such a ultra-light scalar particles are usually postulated to explain dynamic of 
galaxies [12] and are treated as a bosonic cold dark matter candidates [13, 14]. To 
obtain mass of the scalar particle of order 10~ 22 eV one needs 

o; BD ^-^ + 10- 22 , (2.47) 

which is very close to the conformal coupling value [5]. 

Now we show that the ACDM model can also emerge without assumption about 
specific initial conditions but rather as a specific assumptions about the model param- 
eters. For dust matter w m = the equations (2.26) reduce to 

^A/,o=(l-^Wo, (2.48a) 

^n,o = 3 J^~_ 3) ( 4n n ^ - 35 Ax + 8(n - 3) A\\ (-^) ~" , (2.48b) 

n — 4/ \ / an \ -3+n 

^3n,0 = 36{2n -3) ( "" ^ H ~ 3 ) ^^ ~ 35AX ~ 8n AA ) (^) J ' ^- 48C ) 

Expanding these quantities for |n| «C 1 and assuming that n 2 pa and nAx = nAy = 
nA\ pa one obtains 

n„,o ~^ ( $ Ax + 8AA J , (2.49a) 

^Vo -^ (l^i - 35 Ax J (^) ~ 3 , (2.49b) 

and the Hubble's function (2.25) is 

ff HV-i-(n*A«(*rWn«,~A»(.*r > ) (-V 3 . (2.50) 



Note that in the general case of the monotonic approach to the deSitter state we have 

5 4 

-n(n-3) (2.51) 



3 + 2w BD 9 

and the assumption that |n| <C 1 should be interpreted as smallness of the ratio 

5 A , 

n. (2.52) 



3 + 2u; Rn 3 



10 



From (2.22b) we have x(a) ~ Ax (-§y) which gives x(oo) ~ Ax ("%") • From 
the definition of the variable x we obtain that at the present epoch 

<^) = 4j . (2.53) 

H(f> o 

If we assume that in the model we include only the baryonic matter Q m = Qt, m and 
note that in the Brans-Dicke theory the field can be identified as the inverse of the 
gravitational coupling 

<t> = ^ (2-54) 

which is now the function of the spacetime location we obtain the following form of 
the Hubble's function 

Qbm.,0 + Qdm,0 + I ^bm,0 + ^dm,0 I (2.55) 



^H(a )J \ J \ J \a 

where the present dark matter density parameter is 

4 G 



n 



dmfi 



3HG 



o 



(2.56) 



From the Planck satellite data [1] we have that the present value of the density pa- 
rameter of the dark matter Qdm o ~ 0.266 and we obtain 



G 



o 



0.1995. (2.57) 



HG 

Two different assumptions lead to very similar results. 

2.2 Low-energy string theory limit 

The Lagrangian density of the low-energy limit of the bosonic string theory [15-17] 
can be presented in the following form 

£ = e- 2 *( J R + 4V a <l>V a <l>-A) (2.58) 

where $ is the dilaton field. Making the substitution <fi = e -2 * one obtains the Brans- 
Dicke theory with u; BD = — 1 and V(4>) = A<p. Neglecting the matter, the two theories 
are identical, but they differ in their couplings of the scalar field to the other matter 
[18]. 

The constant value of the Brans-Dicke parameter u BD = — 1 leads to 

4 
<J = -n(n-3) (2.59) 

9 



11 



and for dust matter w m = equations (2.26) reduce to 
tt M ,o = ( 1 " 7 12 Q J ^ m ,o , (2.60a) 



n{n — 3) 

fi ^o =-, ^77 ^ ( 3n fi„ M - n(n - 3) Aa; + 6(n - 3) AA ) f-^r) ~" , (2.60b) 

n{n — 3){2n — 3) \ / \a w / 

^,o = ^,^-3) ( " 3(n " 3) fi m>4 - n(n - 3) Aa: - 6n AAJ g) ^ . 

(2.60c) 

In this case, we also expand these quantities for |n| <C 1 which indicates that \S\ ^C 1 
and we are very close to the quadratic potential function. The resulting Hubble's 
function up to linear terms in initial condition is 

§t^t) - 1 ~l ( fim >° - 2x ( fl °) + 4AA ) + l ( Q ™> - 2x ( fl °) + 4AA ) (7) + 

2AA - 4fi m , (-) 'Vf- 
V a o/ / \ a o 

(2.61) 

One notices that in the case of the low-energy string theory limit we obtain a different 
form of the Hubble's function, applying the same type of expansion, but the difference 
lies in the direct connection of the n parameter with the shape of the potential function 
at the critical point. 

The resulting Hubble's function differs by the presence of the term proportional 
to the natural logarithm of the scale factor. At the present epoch this term vanishes 
but in the past it could play critical role. 

3 Conclusions 

Investigating the emergence of the Standard Cosmological Model ACDM from the 
Brans-Dicke theory of gravity shows the significant problem in modern cosmology, 
namely, the problem of initial conditions. Cosmology is not a usual physics of the 
universe because of the presence of this problem at the very beginning. The strategy 
of the modern cosmology is to study all the possible evolutions for all admissible 
initial conditions and then test them against the astronomical data. Which evolutional 
scenario is realised for our universe can be decided by observations and the quality of 
obtained astronomical data is crucial. 

The dynamical systems methods are especially useful in this context. These meth- 
ods give the opportunity to study all the evolutional paths for all admissible initial 
conditions. However cosmology is not just pure mathematics because we need to con- 
strain possible evolutions represented by the trajectories in the phase space by the 
observations, that is we need to find the initial conditions. Cosmology is the science 
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of the evolution of the universe as well as the initial conditions which give rise to this 
evolution. 

In this paper we have shown how the Standard Cosmological Model (ACDM) 
emerges from the Brans-Dicke theory. There are two ways of approaching the ACDM 
model : a monotonic one and an oscillating one. Which specific case takes place cru- 
cially depends on the initial conditions for the universe as well as the model parameters. 
While the leading term in both cases represents ACDM type evolution, there are addi- 
tional terms which describe the asymptotic state of the system (a node or a focus type 
critical point in the phase space). The density parameters for the dark energy as well 
as for the dark matter are emergent parameters from the Brans-Dicke cosmology. The 
values of these parameters are fragile due to value and the sign of the second derivative 
of the scalar field potential at the critical point. 

Finally, from the recent Planck satellite data, we obtain constraint on the vari- 
ability of the gravitational coupling constant as well as the lower limit of the mass of 
the Brans-Dicke scalar field at the deSitter state. 

Our results are independent on the special form of the potential function for the 
Brans-Dicke scalar field and indicate that the value of its second derivative at the 
critical point is crucial for the expansion scenario of current universe. 
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